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I use QCD sum rule ideas to construct models for generalized parton distributions. To this end, the perturba- 
tive parts of QCD sum rules for the pion and nucleon electromagnetic form factors are interpreted in terms 
of GPDs and two models are discussed. One of them takes the double Borel transform at adjusted value of 
the Borel parameter as a model for nonforward parton densities, and another is based on the local duality 
relation. Possible ways of improving these Ansatze are briefly discussed. 



1 Introduction 

The concept of Generalized Parton Distributions Q] [2] [3] |4) is a modern tool to provide a more detailed 
description of the hadron structure. The need for GPDs is dictated by the present-day situation in hadron 
physics, namely: i) The fundamental particles from which the hadrons are built are known: quarks and 
gluons. ii) Quark-gluon interactions are described by QCD whose Lagrangian is also known, in) The 
knowledge of these first principles is not sufficient at the moment, and we still need hints from experiment 
to understand how QCD works, and we need to translate information obtained on the hadron level into the 
language of quark and gluonic fields. 

One can consider projections of combinations of quark and gluonic fields onto hadronic states |P) 
(Q\q a (zi) qp(z2) \P), etc., and interpret them as hadronic wave functions. In principle, solving the 
bound-state equation H\P) = E\P) one should get complete information about the hadronic structure. 
In practice, the equation involving infinite number of Fock components has never been solved. Moreover, 
the wave functions are not directly accessible experimentally. The way out is to use phenomenological 
functions. Well known examples are form factors, usual parton densities, and distribution amplitudes. The 
new functions, Generalized Parton Distributions UlEEl&l (for recent reviews, see |5]|6)), are hybrids of 
these "old" functions which, in their turn, are the limiting cases of the "new" ones. 

2 Generalized parton distributions 

Generalized parton distributions parametrize nonforward matrix elements of lightcone operators. For, 
example, the twist-2 part of the vector operator built of quark fields O^(z) — ■tjj(~z/2)j ti tfj(z/2) in the 
simplest case of a (pseudo)scalar hadron, e.g., pion can be parametrized in two ways. The first one is in 
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terms of the off-forward parton distribution | DEI El 

(P - T - | 0»(z) \P+ r -)=2P» j dxe~^ Pz *>H(x, £, t) , (1) 

(where £ = (rz)/2(Pz) is the skewness of the matrix element and t = r 2 ) or in terms of two double 
distributions (DDs) JT]|3l|3 

r r f 1 /• 1 ~I /J I 

(P_ \0"(z)\P + -)= / / dae- i/3(Pz) - M(rz)/2 {2P^F(/3,a,i)+r' i G(/3,a,i)} . 

(2) 

The variables a;, £ of OFPDs (/?, a of DDs) can be interpreted as momentum fractions: initial and returning 
quarks carry the momenta (x + and (x - (,)P + , (f3P + + (1 + a)r+ /2 and (3P + - (1 - a)r+/2), 

respectively. The functions H (x, £, i), a, t), G(/3, a, t) are related by 

/■l /•1-I/3] 

H(x,£,t)= d(3 da6(x- (3-£a){F(l3,a,t) + £G(f3,a,t)} . (3) 

J-l ^-1 + 1/31 

The resolution of the apparent discrepancy of describing the same object in terms of one or two functions 
is based on the observation that the choice of two DDs F, G is not unambiguous QIH): one can perform 
transformations which do not change the combination dF/d(3 + dG/da |8|. In particular, there exists a 
DD representation in terms of a single function (9) 



d(3 / dae- tl3{Pz) - ia{rz)/2 {2f3P^ + ar^}h{(3, a,t) 
-l J-i+\p\ 

f) fl rl-\P\ 

= 2i^— d/3 dae^ {Pz) ~ la{rz)/2 h{(3 ia ,t) . (4) 



dz^ j_ x j-i+|/3| 

The generalized parton distribution functions provide a very detailed description of the hadronic structure. 
They include information contained in simpler functions, like usual parton densities f(x) and form factors 
F(t), reducing to them in particular limits. The forward limit gives H (x. £ = 0, t = 0) = f(x), while the 
local one produces the reduction formula 



J H(x,£,t)dx = F(t) 



(5) 



Equivalent relations between DDs, parton densities and form factors also can be written. 

Intermediate in complexity are nonforward parton densities 1 10 1 J-{x, t) = H(x, £ = 0, t), or GPDs at 
zero skewness. They reduce to parton densities for zero t, and give form factors after integration over x. 
The functions J 7 (x, t) can be also obtained from the F-DDs by integration over a: 

/•i-M 

!F(x,t) = / F(x,a,t)da. (6) 

J-l+\x\ 



Note, that the a-integral of G-DDs is zero 1 1 1 12 1 because they are odd functions of a. Interplay between 
x and t dependence of T(x, t) is an interesting and nontrivial problem. In particular, it is closely associated 
with the question 1 13 1 of interrelation between large-i behavior of hadron form factors and x — * 1 shape 
of parton densities. 

GPDs accumulate information about long-distance interactions, hence, they are nonperturbative func- 
tions. Possible ways to get theoretical estimates for them include lattice QCD 1 14| and QCD-inspired 
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models 1 131 1161 IT7II18I . Building self-consistent models of GPDs is, however, a rather difficult problem, 
because one needs to satisfy many constraints which should be obeyed by GPDs. They include spectral 
properties, polynomiality condition, positivity, relations to parton densities and form factors QEllll, soft 
pion theorems 1 19 20 1. Most of these conditions are satisfied, of course, in perturbation theory, and there 
were attempts 1 21 22 1 to use perturbative expressions for modeling GPDs. Still, there remains a question 
about relation of perturbative results to functions describing nonperturbative dynamics. In this respect, 
QCD sum rules [23 1 look as an attractive possibility, being an approach which is closely related to Feyn- 
man diagrams. Its basic concept, quark-hadron duality, provides a tool for translating perturbative results 
into statements about nonperturbative functions. In the past, QCD sum rules were used to get information 
about form factors 1241 1251 l26l 1271 f28 1 and parton densities (see |29| and references therein, and |30|). 
A natural idea is to apply the QCD sum rule techniques to model GPDs. This idea in the pion case was 
already elaborated in some detail by the Bochum/Dubna group 1 1 1 1. My goal in the present paper is to give 
interpretation of the approach in terms of double distributions, to discuss the nucleon case, and to combine 
the idea with recent developments. 



3 QCD Sum Rule for Pion Form Factor 

Basic objects of QCD sum rule analysis 1 23 1 are correlators of local currents with quantum numbers of the 
hadrons one intends to study. The usual choice for the pion is the axial current j a = d^j^ a u. Its projection 
on the single-pion state (0|j Q (0)|p) — ifnPa is specified by the pion decay constant f n . To study the pion 
form factor, one should consider correlator of three currents 1 241 1251 

T^(PuP2) = i J e-^+^(0 | T{j p (z a )Ji'(p)j+(z 1 )}\ 0)d 4 Zl d 4 z 2 , (7) 

where J M = e u u^u + edd^^d is the electromagnetic current. The pion-to-pion transition term corre- 
sponds to 

(0\jf,(z 2 )\p 2 }{p 2 \J^0)\Pi)(Pi\^(^)\0) , 
where the pion form factor contribution 

(P2|J"(0)|pi) = 2P»F 7I (Q 2 ) 

appears in the middle matrix element. As usual, Q 2 = —q 2 or Q 2 = —t, if the GPD notation t is used. 
The relevant invariant amplutude can be extracted by taking the projection 1 25 1 

T^p llP2 ) ee n a n^(p llP2 )/(nP) 2 , 

with n a chosen to be a lightlike vector with equal projections on p\ andj>2, (npi) = (np 2 ) = {nP). Since 
(nq) — and n 2 — 0, the projection kills the structures containing q a ,qp and g a p. Still, the projection 
may contain terms, which cannot be directly related to the form factor contribution. In what follows, 
we will always omit them without explicit notice. 

The starting point of the QCD sum approach 1 23 1 is the dispersion relation for invariant amplitudes that 
appear in the correlator, in fact, the double dispersion relation in case of a three-point function I24l l25l . 

T(p 1 ,p 2 ,Q ) = -o / ds i / ds 2 7 2T7 2\ ■ ( 8 ) 

tt 2 7o Ja (si -PV(S2 -P 2 ) 

One should find the expression for the perturbative version of the correlator and nonperturbative corrections 
which modify the spectrum in the p 2 channels, converting the free-quark spectral density into a function 
containing physical hadrons. From a practical point of view, it is more convenient to consider the double 
Borel transform 1241 1251 

-i rOO f>OC 

$(ti,t 2 ,Q 2 ) = ^ / dsi ds 2 p( Sl , S2 ,Q 2 )e- s ^- s ^ (9) 
77 Jo Jo 
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in which power weights are substituted by the exponential ones. Formally, the action of the Borelization 
operator is given by 

B(p 2 ^r){l/(s-p 2 )} = e~ ST . 
QCD sum rule for the pion form factor then has the structure 

f 2 F^(Q 2 )e-< {Tl+T2) + higher states 

Jo Jo 

+ T 2 A( n /T)(a s G 2 ) + t 3 B( Ti /t,Q 2 ) a s (qq) 2 

+ higher condensates , (10) 

where r = n + r%. We kept the pion mass on the left hand side, but it will be neglected from now on. 

4 Structure of perturbative term 

In the case of T M (pi , p 2 ) amplitude, the lowest order perturbative term is given by a triangle diagram. It is 
very convenient to write it in the a-representation (some details can be found in Refs. | 4|l31l ) 



T^ crt (p 1 ,p 2 ) = ^2 / d ai da 2 da. 



1 



0:3(0:1 + o 2 ) 



2P^ 



exp 



a>i — ol-x 

Q 2 a\a 2 +p\a\a z + p\a 2 az 
01+02 + 03 
where A = ct\ + 012 + 0:3. Using the formula 125 1 1261 

B(j>1^T i ){e A ^} = 5{n-A i ) 
one can obtain the double Borel transform 

$"(n,T 2 , Q 2 ) = B( P 2 - n)B{p% -> t 2 )t»(pIpIq 2 ) 

of the perturbative amplitude 

l2 



(ID 



$ Mpcrt (ri,r 2 ,Q 2 ) = £sf~± da ld a 2 da 3 ^+^) + ^oi_o. 



Sin 



0103 
A 

It is instructive to rewrite this expression using new variables 



S(T2- 



a 2 a 3 



A ) 



exp 



-Q 



2 OL\U 2 



A 



(12) 



o 3 oi _ a 2 _ 

x = — ; — = p(l - x) = px ; — = (1 - p)x = px . 



This gives the following integral representation 

^ pert (ri,r 2 ,Q 2 ) = -^ 



zcfe /" dp x(l - x) {2xP fl + x(p - p)^} 



dA 

Jo Jo 

S (n — pAxx) S (t 2 — pAxx) exp [— Appx 2 Q 2 ] 

Two delta functions can be used to perform integration over A and p. The result is 



(13) 



^ pert (ri,r 2 ,Q 2 ) = 



2tt 2 (ti + t 2 ) 7o 



x cfe ^ 2xP M + xg' 



» T i - r 2 
n + r 2 



exp 



_ Q2 (l-*)nT 2 



x(Tl + T 2 ) 



(14) 
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5 GPD interpretation of perturbative term 



In case of GPDs, we should substitute the local current by the bilocal operator O ll (z), in which quark 
fields are separated by a lightlike distance z, and then use parametrization in terms of OFPDs or DDs. In 
fact, nothing prevents us from taking z equal to the projection vector n. Then calculation of the triangle 
diagram contribution completely parallels that for the purely scalar case discussed in Refs. I4ll21l . The 
conversion to DD variables is especially transparent. One should just use the fact that the spinor factors 
/ k 1 , jk 2 corresponding to the numerators of quark propagators adjacent to the composite vertex can be 
written in the a-representation with vectors k% , hi given by 



, a 3 a 2 „a 3 

kl=Pl ~ +r - = P - 



1 



C*2 - a l 

A 



r a 3 a x a 3 



i 



«2 - Ot-l 

A 



r 

2 ' 



while the corresponding quark momenta in DD variables are 0P + (1 + a)r/2 and 0P — (1 — a)r/2, 
respectively. Thus, a 3 / A should be interpreted as 0, and (a 2 — c*i)/A as a. This mnemonics helps to 
understand the DD representation of the triangle diagram 



T^(p uP2 ;a,0) = J J_ da lC la 2 da 3 6 (j3 - y) 



a 3 (ai 



OL2) 



A 2 



2P^ 



A 2 
a 2 — a.\ 
A 



exp 



a 2 — O-i 
A 



-Q 2 a\a 2 +pjaia 3 + p 2 a 2 a 3 
ai + a 2 + a 3 



(15) 



It differs from the representation for T^(j>i,p 2 ) by two delta functions relating the DD variables a, (3 to 
the a-parameters. Note also, that defining GPDs we treat the momentum transfer r = pi — p 2 as "going 
upwards" in the i-channel, i.e., we take r = —q. Using the delta functions to eliminate two integrations, 
we obtain 



T^{pi,p 2 ;a,(3) 



36(0) 
4tt 2 



0(1-0) {20P' 1 + ar lx } 



dXe 4 



j^Q 2 [(!-/3) 2 -a 2 ] 



x exp <^ -A/3[(l -0- a)pj + (1-0 + a)p 2 2 ] 



(16) 



The restriction > reflects the obvious fact that the triangle diagram involves only valence quarks. Now 
it is straightforward to calculate the double Borel transform 



^ pert (ri,r 2 ;a,/3) 



3 0(0) 



2^(ti+t 2 ) 

T2 - Tl 



(1-0){20PI J + ar»} 



(1-0) 



Tl + T 2 

Note, that the spectral property \a\ < 1 — is manifest in this expression. 



(17) 



'Borel" model for nonforward densities 



Integrating DD over a, we should get nonforward parton density (see Eq.(|6j). Performing the integral, we 
obtain 



1-/3 



3 6(0) 



^^ t (Ti,r 2 ;a,0)da = , 
1+13 2tt 2 (ti+t 2 ) 



(l-0)\20P» + (l-0) 



T~2 ~ Tl 
Tl + T 2 



x exp 



(1 - 0)Q 2 TlT 2 
0(Tl + T 2 ) 



(18) 
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Comparing this result with Ea.(ll4>. we see that the integration variable x in that equation has the meaning 
of the momentum fraction, and the integrand can be treated as a nonforwardparton density J 7 (x, t = —Q 2 ). 

Writing QCD sum rule for the pion form factor, we should treat symmetrically both pions, the initial 
and the final one, i.e. take the Borel parameters ti , r 2 equal to each other, as it was done in Refs. 12411251 . 
Such a choice corresponds to DD ^ proportional to 5(a). In other words, each quark takes exactly a half 
of the t-channel momentum r: there is no spread in the distribution of r among constituents. Furthermore, 
the G part of DD vanishes for ri = r 2 . As a result, the perturbative term for the off-forward distribution in 
this case 

HZ eT \x,t,t)\ T1=T ^ r = 6N a x(l - x) exp 

coincides with that for the nonforward parton density, or in the forward t — limit, with the perturbative 
term for the usual parton density f(x). The perturbative term suggests 6x(l — x) for the shape of the 
normalized parton density. This does not look realistic, even for a valence distribution. However, as 
shown in Ref. |30|, nonperturbative corrections shift the maximim of the distribution to smaller x, and 
then DGLAP evolution from a low normalization point /.i 2 ~ 0.25 GeV 2 (to which QCD sum rules refer) 
produces acceptable valence distributions for the pion. In general, it is a rather popular idea that there 
are no skewness effects at a low normalization point, and one can start with the forward approximation, 
generating nontrivial £ dependence through evolution I32i l33l . 

Our result dl9l gives an example of a nontrivial interplay between x and t dependence of a nonfor- 
ward parton density. It has the same form as the result of calculation of the overlap contribution of two 
lightcone wave functions ^(x,k±) with the Gaussian ~ cxp[— a 2 k 2 _/x(l — x)] dependence on trans- 
verse momentum. The parallel between the Borel transform and Gaussian wave functions is well known 
(see, e.g., Refs. 1261 1341 ~). It can be explained by the exponential weight e~ ST and the lightcone form 
s = (fc^ + m q ) 2 /x(l — x) for the invariant mass of the qq system. This analogy with the wave function 
description suggests to take the Borel transform <j>^P ert (r, r; a, j3) at a particular (adjusted) value of the 
Borel parameter as a model for the pion DD, treating r as the width scale of a Gaussian wave function. 
From the QCD sum rule point of view, the model 

-i />oo />oo 

f 2 F*(Q 2 ) = — ds x d S2 pP° rt (si,s 2 ,Q 2 )e- (si+S2)r (20) 
77 Jo Jo 

corresponds to taking such a value of r for which the condensate corrections and the subtraction of contri- 
butions due to higher states perfectly cancel. For each fixed Q 2 , such a value of r exists, but in principle it 
may depend on Q 2 . The absence of such a dependence can be expected only if the description of the pion 
vertex by a Gaussian wave function is a very good approximation of reality. However, inside the QCD sum 
rule approach, the wave function backing of this model faces difficulties. In particular, one can try to check 
the normalization condition -F ff (0) = 1 by using the Ward identity relation 

P pcn ( Sl , s 2 ,Q 2 = 0) = 7rS( Sl - s 2 )p pcrt (si) . 
between three-point and two-point function densities. The resulting expression 

i r°° 

f%F*(0) = - p^\s)e- 2sT ds (21) 

I" Jo 

matches the expression for f 2 derived from the two-point sum rule 

(f*) 2 = ~ p peit (s)e- sr ds (22) 

71" Jo 

only if one decreases the r parameter of the three-point function by factor 2 compared to that used in the 
two-point function relation. 



2x 



It 



(19) 
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7 Local quark-hadron duality model 

Another approach to obtaining predictions for hadronic characteristics from their analogues calculated for 
free-quark systems is suggested by the local quark-hadron duality hypothesis. For the pion form factor, it 
gives 1 25 1 



f, 



lF^ D {Q 2 ) = \ f° dsi r ds 2 /5 pcrt (si,s 2 ,<2 2 ) , (23) 

T Jo JO 

where sq is the duality interval. This relation corresponds to the r = limit of the full QCD sum rule 
l llOt in which the higher states are modeled by the perturbative spectral density starting at sq in both si 
and S2 directions. When r — > 0, the condensate corrections vanish, and the exponential weight e~( sl + S2 ) T 
converts into 1. The local duality approach has no problems with the Q 2 = limit. It gives for f 2 F^ D (0) 
the same expression as the local duality relation 

(0 2 = - Pp pelt (s)ds (24) 



based on two-point sum rule gives for f 2 . Numerically, the result p pcrt (s) = 1/Att |23| fixes so at 
4n 2 f 2 ~ OJGeV 2 This value is also the result of fitting full QCD sum rule, with condensates 
included. 

The duality interval sq has the meaning of the effective threshold for the onset of higher states. Since 
the location of higher states is fixed, sq has good chances to be Q 2 -independent. This hope is supported by 
the fact that the local duality prediction [25 35 1 for the pion form factor 

l + 6s /Q 2 \ a s /ir 



(with a s /ir = 0.1 and sq = An 2 / 2 ) is in perfect agreement with the latest Jefferson Lab measurements 



1361 . The 0(a s ) term in Ea. J25t is the simplest interpolation |35 1 between the Q 2 = value fixed by the 
Ward identity and the large-Q 2 asymptotic behavior F 7T (Q 2 ) — > %na s f 2 /Q 2 [[37 38| due to the one-gluon 
exchange. 

Knowing the double Borel transform one can obtain the spectral density by the inverse transformation. 
It is convenient to write the result in a form similar to the lightcone representation [ 39 1 

P pcrt (s 1 ,s 2 ,Q 2 ) = ^- J xxdx J d 2 K± S( Sl -K 2 ± )S( S2 ~(K ± +q ± ) 2 ) , (26) 

where q± is a two-dimensional vector with q\ = xQ 2 /x. From this representation, it is evident that 
the "Borel" model J20i corresponds to Gaussian wave functions, while the local duality d23i corresponds 
to step-like effective wave functions ~ 9{n\ < so) (note, that k± differs from the usual fcj^ by \fxlt 
rescaling). Thus, to get the form factor through the local duality formula, one needs to calculate the area 
of overlap of two circles having equal radii. This gives the representation 1 39 1 




F^ D {Q 2 ) = — [ ^ d WarccosWf^-Wf^( 1 . - ^- ) } (27) 

Its integrand can be treated as the local duality model for the nonforward parton density T(x, t — —Q 2 ). 
Simple structure of Eq. i26\ allows to write form factor also in the impact parameter representation 

-l 2 

dx , (28) 



F^{Q 2 )^^ [ ^L e ^ ±) f\- 



where Ji(z) is the Bessel function. Again, the integrand gives a model for the nonforward density in the 
impact parameter space. 
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8 Nucleon form factors 

QCD sum rule analysis of the nucleon form factors is based on the study of the 3-point correlator 

T^(PuP2) = J(0 | 7^(22)^(0)^)} | 0)e-^^ +l ^d 4 z 1 d 4 z 2 (29) 
of the electromagnetic current J M and two Ioffe currents 77, 77 1 40 1 with the nucleon quantum numbers, 

v = e abc (u a C-\u b ) lpl5 d c . 

Here, C is the charge conjugation matrix, {a, b, c} refer to quark colors, the absolutely antisymmetric ten- 
sor e abc ensures that the currents are color singlets, and a, (3 are Dirac indices. The amplitude T£„ is the 
sum of various structures: P M P = V^(P), P, ie^ Xpa P\<Zp757o- = A fJ, (P, q), Q 2 7 M , etc. To compare 
contributions of different structures, one should specify the reference frame. A natural choice is the infinite 
momentum frame (IMF), where P M is in the plus direction and P + — > 00, while q^ = q^_ is fixed. Note, 
that neglecting q p compared to P' 1 is exactly what we did in the pion case by taking projection on the light- 
like vector n orthogonal to q. The leading IMF structure is clearly V M (P): it does not contain the "small" 
parameter q. For p\ = p 2 , this structure satisfies the transversality condition q t _ l V fJ '(P) = 0. Another 
structure possessing this property is A' i (P, q). These two structures have the most direct connection with 
the P component of the two-point function defined through the correlator (0 | T{rjf3(z)fj a (0)} | 0). In the 
decomposition of the proton-to-proton transition part of the correlator T£g (pi ,p 2 ), the structure (P) is 
accompanied by the nucleon Pi (Q 2 ) form factor, while the (P,q) structure is accompanied by the form 
factor Gm(Q 2 )- The double Borel transform of the invariant amplitude Ty(pf,p 2 , Q 2 ) related to V^(P) 
structure has the form |41 1 



1 f _ 

$i(ti,t 2 ,(2 2 ) = ■ ttt / dx [3e u x 2 - (2e„ - e d )x 3 ] exp 

(2tt) 4 (ti +T 2 yj 



Q 2 XTlT2 



(30) 



Taking t\ = t 2 = t, one obtains the Borel model for the Pi(Q 2 ) form factor. Similarly, the double Borel 
transform of the Ta(p 2 ,p 2 ,Q 2 ) amplitude |41 1 



®m{ti,T2,Q' z )= jz , 4 , = ; 7^ I dx3e u x 2 exp 

(27r) 4 (ri + T2)' 1 



Q2 xnn 



x(Ti + T 2 ) 



(31) 



gives the Borel model for the magnetic form factor Gm(Q 2 )- Using Gm(Q 2 ) = Fi{Q 2 ) + F-iiQ 2 ) we 
obtain 



1 /-i 

$2(1-1, t 2 , Q 2 ) = ■ -rr / dx (2e u - e d ) x 3 exp 



x{n + r 2 ) 



(32) 



for the double Borel transform of the amplitude related to the F 2 (Q 2 ) form factor. The integrands of 
these representations can be treated (for t\ = t 2 = r) as models for the corresponding nonforward parton 
densities. In the Q 2 = limit, one obtains models for forward parton densities. The use of local duality 
changes the exponential factor exp[— xQ 2 t/2x] into some function of xQ 2 /xso without changing the 
pre-factors, i.e., forward parton densities. Let us discuss main features of the models for forward densities. 
• In case of $1, the forward densities correspond to usual parton densities. The model gives: 



f™ od (x)=4x 2 (3-2x) = 4(l-x) 2 {l + 2x) , (x) = Ax A = 4(1 



x 



Just like in the pion case, the condensate corrections and then DGLAP evolution will shift the distributions 
towards smaller x. One may expect that these effects will modify both distributions in a similar way. So, 
it is interesting to compare relative shapes of the model u and d distributions. The main feature is that 
d distribution has an extra power of (1 — x) compared to u distribution. The extra power of (1 — x) in 
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d distribution is a well known phenomenological observation. For comparison, GRV parametrization for 
normalization point /i 2 ~ 1 Ge V 2 can be rather accurately reproduced by HI 01 

fP hcn (x) = 1.89aT°- 4 (l - x) 3 - 5 (l + 6s) , f d (x) phcn = 0.54aT - 6 (l - .t) 4 - 2 (1 + 8a;) , 

with d distribution having extra power (1 — x) ' 7 for x — ► 1. It is easy to check that the ratios 
f™ od (x)/f™ od (x) and f% hen {x)/fv hen (x) are very close to each other. 

• In case of $2, the densities e a (x) correspond to the forward limit of the spin-flip GPDs 

e a (x) =E a (x,t = 0,t = 0) . 

They are inaccessible in deep inelastic scattering and other inclusive processes. However, the normalization 
integrals n a for these functions are related to the anomalous magnetic moments of the nucleons: 



and, hence, they are known: 



— H~ K n — 1.65 , — *2iK n H~ Kp — 2.07 . 

The model densities are 

e™ d (x)=8(l-x) 3 , er d (x) = -4(1 - x) 3 . 

The normalization integrals for these functions are K™ od = 2 and K d nod = — 1. The correct nontrivial sign 
of Kd is a very encouraging indication that the model is a reliable starting point. Another feature of the 
model e Ui d(x) distributions is that they both have an extra power of (1 — x) compared to f u {x). Asaresult, 
the F2(Q 2 ) form factor, as we will see later, decreases faster with Q 2 than Fi(Q 2 ) and Gm(Q 2 ), again in 
agreement with experimental observations. 



9 Form factors at large momentum transfer 

In the Borel model, the form factors are given by the integrals 



F a (Q 2 ) = / dxf*(x)exp 



2.r 



(33) 



involving the model parton density f B (x) and the exponential factor 



-{1-x)Q 2 t /2x 



. Hence, at large Q 2 



t/Q 2 . For the 
< x < 1, i.e., 



the form factors are dominated by integration over the region where Q 2 tx ~ 1 or 1 — x 
local duality model, the whole integration region over x is restricted to (1 + 4so/Q 2 )~ 
1 - x < 4s /Q 2 for large Q 2 . 

Hence, the result of integration in both models is completely determined by the behavior of parton 
densities at x close to 1. Namely, if f(x) ~ (1 — x) v for x — > 1, then the relevant form factor drops like 
l/(Q 2 ) i/+1 at large Q 2 . This gives l/Q 4 forthe asymptotics of the pion form factor, l/Q 6 forthe large-Q 2 
behavior of the nucleon form factors Fi(Q 2 ) and Gm(Q 2 ), and 1/Q S for ^(Q 2 )- All these results seem 
to be in contradiction with the experimentally established exponents of the power-law behavior of these 
form factors, so one may be tempted to conclude that our models have no chance to describe the data. 

In fact, as already mentioned, the local duality prediction (I25> is in excellent agreement with the results 
of recent JLab data. In the nucleon case, the local duality calculation of G P M {Q 2 ) performed in Ref. 1 26l l41l 
agrees with the data up to Q 2 ~ 20GeV 2 (see also Ref. |42|, where the curve based on Refs.|26 41 1 is 
compared with the results of other approaches). The ratio F^ (Q 2 )/Ff (Q 2 ) as calculated in the local 
duality model agrees with the data based on Rosenbluth separation |43|, though at the highest Q 2 it is 
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somewhat lower than the results of the polarization transfer experiments |44|. In all cases, deviations do 
not exceed 20-30%. This means that the model curves mimic the "canonic" Q 2 behavior of these form 
factors (1/Q 2 for F^Q 2 ), 1/Q 4 for if (Q 2 ) and G P M (Q 2 ), and 1/Q 6 for Ff(Q 2 )). 

The resolution of the paradox is based on a trivial observation that the model curves are more com- 
plicated functions than just pure powers of 1/Q 2 . Their nominal large-Q 2 asymptotics is achieved only 
at very large values of Q 2 , well beyond the accessible region. Thus, conclusions made on the basis of 
asymptotic relations might be of little importance in practice: a curve with "wrong" large-Q 2 behaviour 
might be quite successful phenomenologically in a rather wide range of Q 2 . 



10 Improved Ansatze 

The models discussed above, of course, have some drawbacks. In particular, the small- a; behavior of the 
model parton densities is unsatisfactory: it does not have the standard Regge-type behavior J(x)\ x ^q ~ 
x -a(0) jjjg excuse j s suc h a behavior cannot result from a single lowest-order triangle diagram: to get 
it, one should add an infinite number of diagrams with the ladder structure in the Q 2 channel, and perform 
summation of all contributions. A simple solution of this hopeless problem is to take experimental forward 
distributions in Eq. J33i instead of the model ones. Such an approach was successfully used in Ref. ll II 
for the pion form factor and in Refs. 1 1 01 1171 for the proton F[(Q 2 ) form factor. The modified Borel (or 
Gaussian) model for Ff(Q 2 ) was able to successfully describe the data up to Q 2 ~ lOGeV 2 . 

Next question is about the exponential factor. The Regge picture suggests x^ a ^> behavior at small x or 

T{x,t)= f{x)x- {a(t) ~ a(Q)) (34) 
model for the nonforward densities. Assuming a linear Regge trajectory with the slope a', one gets 

F R (x,t) = f(x)x- a,t = f(x)e- a ' tlnx . (35) 
This Ansatz was already discussed in Refs. 1 5 1 and 1 45 1 . It also provides finite mean squared radii 

(r 2 ) = -6a; [ dxf(x)lnx, (36) 
Jo 

which are in good agreement with experimental values |45 1. For large t, experimental data support Drell- 
Yan (DY) relation 1 13 1 and Bloom-Gilman duality |46|. According to DY, if the parton density behaves 
like (1 — x) v , then the relevant form factor should decrease as l/t^ +1 ^ 2 for large t. The simplest idea is 
to attach an extra (1 — x) factor in the exponential |47 1, i.e. to take the model 

Tr° d M = Ux)x- a '^-^ . (37) 

To calculate F%, we need an Ansatz for the spin-flip nonforward parton densities £ a (x, t). One can assume 
the same model 8 Rtn ° (x, t) = e a (x) a j-«i( 1 - a ')* a s for T a (x, t), with possibly a slightly different slope 
a' 2 - To model the forward magnetic densities e a (x), we can use the lesson from the triangle diagram 
calculation that e u (x) has an extra power of (1 — a;) compared to f u (x). In fact, one can take the extra 
factors in the form (1 — x) Va with r/ a 's being fitting parameters. Within this approach, it is possible to get 
a rather good description of all four nucleon form factors |48 1 (see Ref. |49| for a similar analysis). 



11 Discussion 

In this paper, we discussed basics of an approach that uses QCD sum rule ideas to build models for general- 
ized parton distributions. The underlying idea is to consider three-point functions in which the hadrons are 
represented by local currents with appropriate quantum numbers. The necessary nonlocality bringing in 
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parameters having the meaning of the hadron size can be introduced in several ways: by nonzero virtuali- 
ties p 2 of the momenta associated with the currents, by taking the double Borel transformation fromp 2 's to 
Tj's, or through the local duality prescription within the duality square sq x sq. In terms of the basic func- 
tion, the spectral density p(s\, S2 , t), these possibilities correspond to integration with different weights: a) 
(s\ — p\)~ 1 (s2 — p\)~ 2 (producing the original amplitude T{p\ ,p§!*))> exp[— s\T\ — S2T2] (producing 
the double Borel transform $(rx,T2, Q 2 )), and c) 6(si < sq)9(s2 < sq) (local duality prescription). We 
have considered only the lowest approximation for p(si, S2,t) (0(a s ) terms for pion form factor were 
discussed in refs.|28 50 1). One of the expected effects of higher-order corrections is the emergence of 
the Regge-type behavior at small x. Since there is no doubt about this outcome, such a behavior can be 
introduced in the model expressions by using experimental parton densities instead of those generated by 
the lowest order term. 

In higher order diagrams, one would also obtain contributions corresponding to the leading large-Q 2 
asymptotics of perturbative QCD, like the one-gluon-exchange term for the pion form factor. All such 
higher-order contributions are suppressed by a s /ir ~ 1/10 factor per each extra loop. Such a suppression 
is manifest in the local duality result d25l > for the pion form factor. Note, that the local duality prediction is 
in perfect agreement with the data despite the fact that the 0(a s ) term containing the hard gluon exchange 
is insignificant compared to the lowest order term. For the nucleon form factors, the leading pQCD two- 
gluon exchange term has a priori suppression by a factor of 100, so it is unlikely to be relevant at any 
accessible momentum transfer. 

Another pQCD prediction is about the x — > 1 behavior of parton densities. In the nucleon case, the 
leading (1 — x) 3 term corresponds to four-gluon-exchange diagrams and ~ 10~ 4 suppression compared to 
the O(a s ) term. The latter has (1 — x) 2 behavior for u quarks and (1 — a;) 3 for d quarks. Nonperturbative 
effects and DGLAP evolution are undoubtedly capable to shift these densities towards the experimentally 
observed shapes. Furthermore, there is no need to make extra efforts to bring in the relative (1 — x) 
suppression of the d density: it is present in the starting approximation. 

The lowest order term also implies a faster fall-off of the ^(Q 2 ) form factor compared to Fi(Q 2 ). This 
effect results from the extra (1 — x) power of e™ d (x) compared to f™ oA (x). It should be emphasized 
that, in general, the large-Q 2 behavior of form factors in the C(a° ) approximation is completely governed 
by Feynman mechanism, i.e., by x ~ 1 integration, so that the large-Q 2 behavior of form factors is always 
determined by the x — > 1 behavior of the parton densities. The specific correlation pattern between the v 
power in (1 — x) v and n in 1/ (Q 2 ) n depends on the structure of a factor like exp[— Q 2 t(1 — x)/2x)] in the 
Borelmodel. As we discussed, it should be modified to exp[Q 2 a'(l — x) Inx] to impose the n = {v+l)/2 
correlation dictated by the Drell-Yan relation. In contrast, in pQCD the large-Q 2 behavior of form factors 
is governed by configurations in which all the valence quarks carry finite momentum fractions Xi, i.e., the 
l/Q 4 behavior of i^i(<5 2 ) is not a consequence of the (1 — x) 3 behavior of the parton densities. The relation 
n = [y + l)/2 in pQCD is just an accidental correlation between two parameters. In other words, there 
is a correlation between v and n in pQCD because both are determined by the same hard gluon exchange 
mechanism, but there is no causal connection between these two numbers. As noted in the pioneering 
paper [38 1, there is no Drell-Yan/Feynman mechanism in pQCD. 

Summarizing, the gross features of generalized parton distributions are dominated by nonperturbative 
dynamics, and, hence, we need nonperturbative approaches to build models for GPDs. The models mo- 
tivated by QCD sum rule ideas have already made several successful predictions, and they also have the 
advantage of being closely related to perturbative calculations, which allows to satisfy nontrivial con- 
straints imposed on GPDs. This makes the QCD sum rule based approach an attractive possibility for 
building realistic models of GPDs. 
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